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Heavy-light decay constants from the step scaling method* 

Filippo Palombi a 

a "E. Fermi" Research Center, c/o Compendia- Viminale - pal. F, F00184 Rome, Italy 

We discuss results for the heavy-light decay constants in the continuum limit of quenched lattice QCD from 
finite size scaling techniques. We disentangle the simultaneous presence of the different energy scales characterizing 
heavy-light physics by first performing simulations at the unphysical volume Lo = 0.4 fm, and then evolving the 
results towards the infinite volume. We find f B , = 192(6) (4) MeV and f Ds = 240(5) (5) MeV. The approach has 
been developed by the APE group at the University of Rome "Tor Vergata" . 



1. Introduction 

The decay constants fhi of the heavy-light 
pseudoscalar mesons depend upon two different 
energy scales. A first energy scale is the one as- 
sociated with the heavy quark mass rrih, and a 
second one is related to the light quark mass mi, 
or Aqcd- As a consequence, numerical simula- 
tions of heavy-light decay constants require tiny 
lattice spacings, in order to properly describe the 
dynamics of the heavy quark (highly localized), 
and large physical volumes, to accommodate the 
light quark (widely spread). In practice, a direct 
reliable simulation of fhi should demand a lat- 
tice size of around L/a ~ 0(100), which is hardly 
affordable nowadays in quenched simulations as 
well as in future unquenched ones. Different ap- 
proaches are available on the market to face this 
problem [3]. I describe a new method which has 
been developed by the APE group at the Univer- 
sity of Rome "Tor Vergata" |1I2| , called the step 
scaling method at Lattice 2002 conference 0]. 



2. Step scaling method 

The step scaling method is a two step proce- 
dure where the heavy-light decay constants fhi 
are first computed on a small unphysical volume 
with Lq = 0.4 fm, and then evolved towards 
larger volumes in order to remove finite size ef- 
fects. The evolution is simply provided by the 



identity 



fhe(Loo) — fhi{Lo) 



fhe(Lo) fhe(Li) 



(1) 



where Lq < L\ < L 2 ■ ■ ■ The basic ingredient of 
the procedure is the computation of the step scal- 
ing function 



f he (mi,m h ,L k ) 



fhe(me,mh,L k - 



(2) 



L k =sL k _ x 



"talk given at Lattice 2003 



i.e. the ratio of the decay constants at two dif- 
ferent volumes and same quark masses. We fix 
s = 2 as the ratio of subsequent sizes and, af- 
ter two evolution steps, we end up on L 2 = 1.6 
fm, which we take as the "infinite" volume, i.e. 
free of finite size effects. At fixed lattice spac- 
ing, both fhe(Lo) and the c's are simulated for 
a set of different quark masses, measured in a 
RGI quark mass scheme. Then, they are sepa- 
rately extrapolated to the continuum limit. Units 
are fixed through the rg scale |5I6| . The main 
advantage of the method is that fhe(Lo) can be 
simulated at physical heavy quark masses with- 
out big lattice artifacts while the cr's preserve in 
all steps a soft dependence upon the heavy quark 
mass, due to cancellations among volume depend- 
ing terms. This makes the extrapolation to the 
physical values of m^ GI and m^ GI numerically 
safe. Moreover, such values can be obtained by 
a self-consistent application of the step scaling 
method to the HL-meson masses |7I8| . 
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3. Computational framework 

The calculation is set up on the Schrodinger 
Functional scheme, with topological parameters 



T = 2L, C = C" = 0, 







(3) 



C and C being the gauge fields on time bound- 
aries, and 8 being a phase which affects the pe- 
riodicity of the fermion spatial boundary condi- 
tions. We simulate the basic SF correlation func- 
tions, i.e. the axial, pseudoscalar and boundary 
normalization bilinears 



a 6 ^ 



(4) 



(5) 



Fl = -?f? E (0(y)7 5 C i (z)CKy')75Cj( z ')) (6) 



3L 6 

yzy'z' 

In terms of these, a finite size definition of the 
decay constants is given by taking the axial cor- 
relation function in the middle of the lattice 

Z A (l + b A am) F ^ 2) (7) 



fht = 



y/L*M A (T/2) 



and the normalization is provided in terms of the 
effective meson mass, still measured in a finite 
size fashion 



aM A (T/2) 



1 



■In 



F{{T/2- 



FUT/2- 



(8) 



Quark flavors are defined according to their RGI 
masses (scale/finite-size independent), defined as 
follows. We start from the bare quark masses 



1 



in 



wi 



(l/2)(d + dl)F A + ac A d^d F F 



2F P 



(9) 
(10) 



T/2 



and connect them to the RGI quark masses 
through proper renormalisation factors known in 
literature |9TTU| 



m^ GI = Z M 



\ + {b A -b P )- 



am,i 



wi 



(11) 



,R,GI 



= Z M Z[l + b m arrii]r 



z= z^Zp (12) 
z A 



In this way, we get a plethora of definitions, all 
with the same continuum limit and differing in 
0(a 2 ) terms: 



1. m 



RGI 



m 



RGI 



-.RGI _ 



,.RGJ 



2m 



RGI 



n RGI 
l 33 ■ 



3 + i 



After choosing a target pair (m 



RGI 



,RGI 



tunc the hopping parameters (ni,Kj) such that 
the values of the masses coming from the vari- 
ous definitions equate the chosen pair. Of course, 
each choice of Kj) produce a different value of 
fht, all the values differing in 0(a 2 ) terms, as it 
can be seen in Fig.l. 



xf 01 /dof = 0.89 




















1 






• 


i - 






' n 






i---f~ » — 


















i 










i 





.71 : - * 
33 - * 

.37 ^ -B- 

.37 

.33 i • - 

.37 hA- - 



0.001 0.002 0.003 0.004 0.005 0.006 

Figure 1. Continuum extrapolation on the small vol- 
ume of /b 3 (Lo). Units in GeV. 

4. Results at L = 0.4 fm 

Simulations of the decay constants on the 
smallest volume Lq = 0.4 fm have been performed 
at five different lattice spacings, using the ge- 
ometries 24 x 12 3 , 32 x 16 3 , 40 x 20 3 , 48 x 24 3 
and 64 x 32 3 . For each lattice spacing a set of 
eight RGI quark masses have been simulated, two 
around mf G/ , two around mf G/ , one additional 
at 4.0 GeV, and three light quark masses at 140, 
100 and 60 MeV. We fix the light quark mass at 
the physical value of mf GI by interpolation and 
we get 



f B ,(L ) = 475(2)MeV 



(13) 



3 



/D s (io) = 644(3)MeV 



(14) 



The errors quoted at this stage are statistical 
only, and they have been evaluated by a jack- 
knife procedure. The continuum extrapolation 
has been obtained by performing a combined lin- 
ear fit in (a/ro) 2 over the last three sets of points 
of Fig.l. We estimate a systematic error of 1%, 
related to this extrapolation, by the discrepancy 
observed after restricting the fit to the two sets 
of points nearest to the continuum. The results 
obtained are clearly unphysical, as they have to 
be corrected through the evolution factors given 
by the cr's. 

5. Evolution steps 

The first evolution step connects Lo to L\ = 
2I/o = 0.8 fm. Three different lattice spacings 
have been simulated in order to perform the con- 
tinuum extrapolation, corresponding to the ge- 
ometries with Lq/ci = 8, 12, 16 and consequently 
Li/a = 16,24,32. The values of the simulated 
quark masses have been halved in order to leave 
the lattice artifacts unchanged. In Fig. 2 the con- 
tinuum step scaling function a(Lo) is plotted vs. 
the inverse of the heavy quark mass with the light 
quark mass fixed at mf GI . At this stage we get 
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Figure 2. Continuum step scaling function a(Lo) as 
a function of l/m^ GI at mf GI = mf GI . Units in 



GeV. 

o-b.OLo) =0.417(3) 
o-u.(L ) =0.414(3) 



(15) 
(16) 



Analogously, the second evolution step connects 
L\ to L2 — 2Li — 1.6 fm. Also here, three lat- 
tice spacings have been simulated with the same 
lattice geometries as before. Again the values 
of the simulated quark masses have been halved. 
Fig. 3 shows the continuum step scaling function 



cr(Li) vs. the inverse of the heavy quark mass at 
rgi _ 
1 



m 



mf GI . We find 
a Bs (Li) = 0.97(3) 
a Ds {Li) =0.90(2) 



(17) 
(18) 
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Figure 3. Continuum step scaling function ct(Li) as 
a function of l/m^ GI at mf GI = mf GI . Units in 
GeV. 

6. Results 

Combining the results of the small volume Lq 
with the evolution factors, we obtain 

f Bs = 192(6) (4) MeV (19) 

f Ds = 240(5) (5) MeV (20) 

The first error is statistical while the second one 
is systematic and due to the uncertainties on the 
continuum extrapolation, on the scale and on the 
renormalisation factors. 

REFERENCES 



M. Guagnelli et al., Phys. Lett. B546 (2002) 237, 
hep-lat/0206023 

G. M. de Divitiis, M. Guagnelli, 

F. Palombi, R. Petronzio and N. Tantalo, 

arXiv:hep-lat /0307005 

Andreas S. Kronfeld, these proceedings. 

N. Yamada, arXiv:hep-lat/0210035 

S. Necco and R. Sommer, Nucl. Phys. B 622 

(2002) 328 



10. 



M. Guagnelli, R. Petronzio and N. Tantalo, Phys. 

Lett. B 548 (2002) 58 

N. Tantalo, these proceedings. 

G. M. de Divitiis, M. Guagnelli, 

R. Petronzio, N. Tantalo and F. Palombi, 

arXiv:hep-lat/0305018 

S. Capitani, M. Guagnelli, M. Luescher, S. Sint, 
R. Sommer, P. Weisz and H. Wittig, Nucl. Phys. 
Proc. Suppl. 63 (1998) 153 
J. Rolf and S. Sint [ALPHA Collaboration], 
JHEP 0212 (2002) 007 



